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As a simplified model of randomly pinned 
vortex lattices or charge-density waves, we 
study the random-field XY model on square 
(d = 2) and simple cubic (d = 3) lattices. We 
argue, and confirm in simulations, that the 
spacing between topological defects (vortices) 
diverges more strongly than the Imry-Ma 
pinning length as the random field strength, 
H } is reduced. For d = 3 the data are consis- 
tent with a topological phase transition at a 
nonzero H c to a vortex-free pinned phase. 
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There is considerable current interest in 
the effects of quenched disorder on ordered 
phases with continuous symmetry. Exam- 
ples are vortex lattices in type-II supercon- 
ductors [1], spin- and charge-density wave 
systems [2] subject to random pinning, as 
well as amorphous ferromagnets with random 
anisotropy [3] and liquid crystals in porous 
media [4]. The random pinning induces con- 
tinuous, elastic distortions of the ordered 
state and it may also induce plastic deforma- 
tion due to topological defects, such as dis- 
locations, which do not represent continuous 
distortions of the ideal ordered state. The 
distortions induced by random pinning in the 
absence of topological defects can be treated 
within elasticity theory and have received 
some attention over the years [5-9]. Here we 
focus instead on the production of topological 
defects by the random pinning [10]. 

The simplest system in which to study 
these issues appears to be the ferromagnetic 
XY model with a random field. Here the 
long-range order in the pure system is ferro- 
magnetism, the pinning is due to the random 
field, and the topological defects are vortices 
in the magnetization pattern. The hamilto- 
nian we consider is 
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where the first sum is over all nearest- 
neighbor pairs of lattice sites, S; is a unit- 
length, two-component (XY) spin at site z, 
and the static random fields h; have rms mag- 
nitude H. The ground state of (1) evolves 
from ferromagnetic and vortex-free for H = 
to a state with all spins aligned with the ran- 
dom fields, and thus a dense array of vortices, 
for large H . In the context of vortex lattices, 
the H = limit corresponds to an unpinned 
Abrikosov lattice, while large H corresponds 
to a vortex-glass ground state at strong pin- 
ning [1,11]. 

Let us first consider small random field H } 
following Imry and Ma [6] . Treating the ran- 
dom field perturbation, the long-range 
ordered state is stable at small H only for 
dimension d > 4. For d < 4 the static elas- 
tic strains in the ground state induced by 
the random field are of order one at the pin- 
ning length £p ~ if~ 2 /( 4 ~ d ). The behavior at 
longer scales is not accessible by simply per- 
turbing in H . This perturbative treatment 
only considers continuous elastic distortions 
of the uniform state, not vortices, which are 
nonperturbative. Let us ask about the sys- 
tem's stability at small H and length scale 
L < ip to static vortex/antivortex pairs in 
d = 2 or to vortex loops in d = 3. The vor- 
tices permit the system to align better with 
the random field on length scale L. This, 
naively, lowers the random field pinning en- 
ergy by order HmL d l 2 , where m is the mag- 
netization density at H = 0. However, the 
added energy of the elastic strains around the 
vortices is order KL^ d ~ 2 ^ \ log(iv)|, where K is 
the spin stiffness. Even at L = £p, the elastic 
energy cost is larger than the pinning energy 
gain by a factor of | \og(H)\. Hence, the sys- 
tem appears to be stable against vortices for 
small H at length scales less than or of order 
^p(H). Thus, we conclude that if the system 
has vortices at small H } the length scale at 
which the vortices appear, <^y, is larger than 
the pinning length, £p. Other than this lower 
bound, how £y depends on H and when and 



how the static vortices proliferate at scales 
L > £p is not known. This is what we inves- 
tigate in this paper. 

If vortices are forbidden, then the long- 
distance spin-spin correlation function is ex- 
pected to decay as a power-law [9]. This 
behavior should also apply in the distance 
range between £p and £y where the system 
is strongly pinned but still largely vortex- 
free. Thus the true correlation length, be- 
yond which the correlations decay exponen- 
tially, should be of order <^y. There are two 
possible behaviors for ^v(H): (i) it may di- 
verge only at H = 0, but with a stronger 
divergence than £p, or (ii) it could diverge 
at a nonzero critical field, H c . Our results 
below are more consistent with the former 
for d = 2 and the latter for d = 3. In the 
latter case there is an intermediate pinned 
phase for < H < H c that is vortex-free at 
the largest length scales, and therefore has 
topological long-range order. This novel or- 
dered pinned phase, if it exists, has power-law 
long-distance spin-spin correlations distances 
and is separated from the disordered, plas- 
tic phase at higher H by a topological phase 
transition at H c where large-scale static vor- 
tices first appear. 

We have performed Monte Carlo simula- 
tions of the random-field XY model (1) on 
simple cubic (d = 3) and square (d = 2) lat- 
tices. In both cases we find that, as expected 
by the above argument, the spacing between 
vortices diverges more strongly with decreas- 
ing H than the pinning length £p obtained 
from the Imry-Ma estimate. For d = 3, the 
dat consistent with a transition to a 

vortex-free phase at a nonzero critical field 
H c > 0. For d = 2, on the other hand, the 
vortex density is better fit by a power-law 
in H } indicating that there is no intermedi- 
ate vortex-free phase. For d = 3 we have 
tried to more precisely locate the proposed 
phase transition at H c using various forms of 
finite-size scaling involving moments of the 
magnetization distribution and the derivative 
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of the magnetization with respect to H } but 
have been unsuccessful in obtaining convinc- 
ing one-parameter scaling {L/£(H)}. This 
lack of good finite-size scaling may be because 
there are two distinct diverging length scales 
in this problem, namely £p and £y, which 
complicates the finite-size scaling. 

Here we report on simulations of large 
lattices (up to 10 6 spins) in the field range 
that we could equilibrate. We simulated 
two copies (replicas) of each sample, one 
with a ferromagnetic initial condition and 
one with the spins initially aligned with the 
random field at each site. We only use 
late-time results where both replicas give 
the same time-independent averages. At 
the lowest fields studied, this required 10 5 
Monte Carlo steps per spin. We studied 
the model with Gaussian-distributed ran- 
dom fields [hi] = 0, [hi • hj] = H 2 6 tr We 
measured the time-averaged magnetization, 
m; =< S; >, at each site. The angle between 
the vectors m; on each nearest-neighbor pair 
of sites was obtained, with the convention 
that it lies between —tt and tt. For each ele- 
mentary square plaquette, these angles were 
added to obtain the total rotation of the 
magnetization on moving around the plaque- 
tte. This sum is a multiple of 2tt; if it is 
nonzero, then there is a vortex in that pla- 
quette in the equilibrium (static) magneti- 
zation pattern. We also measured the cor- 
relation function g(r) = [m; • mj], averaged 
on pairs of sites at distance r along a lattice 
axis. We want to study the ordered-phase 
(low-temperature) behavior, but be at a high 
enough temperature that we can equilibrate 
in not too much computer time. For d = 3, 
where the critical point in the absence of the 
random field (H = 0) is at T C 3D = 2.2 [12], 
we examined T = 1.5, which is enough be- 
low T C 3D to avoid the H = critical regime. 
For d = 2, where T 2D = 0.9 [13], we studied 
T = 0.7. 

The fraction, _/y, of elementary square 
plaquettes occupied by vortices is shown vs. 



H in Fig. 1. The solid and dashed lines indi- 
cate what the slopes would be if the intervor- 
tex spacing was the pinning length, £p, given 
by the Imry-Ma argument, so fy ~ ^p 2 . The 
above argument says that fy should vanish 
more rapidly than this with decreasing H } 
and the data support this. For d = 3, the 
vortex density is not well approximated by a 
power of H over any substantial field range. 

In Fig. 2a we show the correlation func- 
tion for d = 3, which is very well fit by a 
simple exponential: g(r) ~ exp(— r/£). The 
measured correlation length, £, also diverges 
faster than the Imry-Ma £p as H is decreased. 
We find that the vortex density, _/y, equili- 
brates much earlier in the simulation than 
the correlation function, hence we have re- 
sults that we trust for fy to lower field than 
for g(r). 

One naively expects that if there is a tran- 
sition to a vortex-free phase for d = 3 at a 
nonzero H c (which depends on T), the vor- 
tex density fy vanishes and the correlation 
length £ diverges as a power of (H — H c ). In 
Fig. 3 we show that our data for T = 1.5 
are consistent with such a critical behavior 
with H c = 1.35. The apparent exponents 
with H c = 1.35, indicated by the solid lines 
in Fig. 3, are £ ~ (H - H c )~ u with v = 0.85, 
and f v ~ (H - H c ) p with p = 1.4. Note 
that the range of the scaling fits for all the 
quantities in Fig. 3 is less than one decade, 
so this apparent scaling should not be taken 
too seriously. But we can definitely say that 
the data for fy in this field range are very 
different from a power-law critical point with 
H c = and, with the data for £, consistent 
with a power-law critical point with H c near 
1.3. 

For d = 2 the Imry-Ma argument gives 
~ ^/H at T = 0. At finite tempera- 
ture, the magnetization at scale L, m(L), is 
renormalized by thermal fluctuations in the 
critical phase below the Kosterlitz-Thouless 
temperature Tkt- When the Imry-Ma argu- 
ment is modified to take this into account 
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one obtains £p ~ i7~ 2 /( 2- ^). If we then take 
the largest value rj = 1/4 at the Kosterlitz- 
Thouless transition temperature [14], we ob- 
tain the slope 16/7 indicated by the solid line 
near the d = 2 data in Fig. 1. As expected, 
the vortex density for d = 2 is found to vanish 
even faster than this with decreasing H. This 
again indicates that, as found in d = 3, the 
vortex spacing diverges more strongly than 
the pinning length given by the Imry-Ma ar- 
gument. However, for d = 2 the low field 
behavior is quite consistent with a power-law 
in H, just with a larger exponent (roughly 
fv ~ H 3 ) } rather than a phase transition 
at nonzero H c . The correlation function for 
d = 2, (Fig. 2b) and low field does not fit 
a simple exponential as found in d = 3, so 
it is not obvious how one should define the 
correlation length. 

How should one think about the physics 
at length scales beyond £p? One possibility is 
the following: Consider d = 2, for concrete- 
ness. First, forbid vortices and find the lowest 
energy state. This state is pinned with some 
particular nontrivial spin pattern. Now con- 
sider a patch of linear size L, ( L 2 ) with 
L ^> £p. Introduce a vortex- antivortex pair 
in this patch and choose their locations and 
the spins' orientations to minimize the en- 
ergy with this pair present. For L < this 
new energy is presumably typically higher 
than the lowest-energy vortex-free state, but 
for L > it is typically lower, so the true 
ground state has typical vortex spacing £y. 

What does the ground state with vortex 
separation £y look like? The relative spin 
orientation between it and the lowest-energy 
vortex-free state must rotate by 2tt on encir- 
cling any vortex. But the system (for £y ^> 
£p) is strongly pinned, so it will typically cost 
a lot of energy to locally rotate the spins away 
from their local ground state. Thus we ex- 
pect that the relative phase difference will 
typically be concentrated in line defects (like 
Sine-Gordon solitons or domain walls) that 
extend from vortex to antivortex. Thus to 



find the ground state one must optimize not 
only over the positions of the vortices but also 
of these line defects (these defects are present 
only in the relative phases and are presum- 
ably of width of order £p). Since the line de- 
fect is not permitted in the vortex-free state, 
it can have negative energy relative to the 
lowest-energy vortex-free state once its posi- 
tion is optimized. This can then "pay for" the 
positive energy of the vortex cores. This pic- 
ture provides an energy-balance mechanism 
for a phase transition driven by a prolifera- 
tion of topological defects. The description 
of the physics at scale L > £p in this picture 
relies heavily on the existence of topological 
defects (vortices). It is not at all a renormal- 
ized elastic medium of the sort considered by 
Giamarchi and Le Doussal [9]. It is still un- 
clear how to This scenario may also apply to 
random- anisotrop y XY models [3,15]. 

We now comment on the connection be- 
tween the d = 3 random field XY model 
and the Abrikosov vortex lattice with uncor- 
rected random pinning. 

There are two aspects of the vortex lattice 
that are not captured by the XY model. The 
single continuous degree of freedom in the XY 
model, the spin orientation, plays the role of 
the vector displacement of the vortex lattice. 
One consequence of this simplification is that 
the topological charge of a vortex in the XY 
model is a scalar (how much the spin orienta- 
tion winds upon encircling the vortex), while 
that of a dislocation in the vortex lattice 
is the two-component Burgers' vector. The 
other simplification is that the superconduc- 
tor has the additional U(l) symmetry of the 
complex scalar Ginzburg-Landau parameter 
if) under rotations in the complex plane; this 
symmetry is absent in the random-field XY 
model. This U(l) symmetry can be sponta- 
neously broken in the absence of any vortex- 
lattice order in the superconductor, yielding 
the vortex glass phase [1,11]. The random- 
field XY model has no analogous ordered 
phase. Another consequence of the U(l) sym- 
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metry is that in the vortex lattice an inter- 
stitial (or vacancy), being a change in the 
vortex number, is itself a topological defect; 
the random-field XY model has no analogous 
defect . 

The vortex-free pinned phase at interme- 
diate disorder H that we discuss in this paper 
for the random field XY model corresponds to 
a pinned, superconducting, dislocation-free 
Abrikosov vortex lattice phase in a type-II 
superconductor [11]. In the latter system, 
the structure factor, S'(q), would have power 
law singularities at the basic reciprocal lattice 
vectors, Q, of the form S(q) ~ |q — Q|~( 2_?? ) 
[9]. By increasing the disorder H in the 
random-field XY model, we drive the sys- 
tem into the fully disordered phase. In that 
case, the disordered phase is not thermody- 
namically distinct from the high-temperature 
paramagnetic phase. For the superconduc- 
tor, on the other hand, there are at least two 
noncrystalline phases: the superconducting 
vortex glass phase at low T and the resistive 
vortex liquid phase at higher T. The transi- 
tion directly from the vortex lattice phase to 
the vortex liquid involves loss of both crys- 
talline and superconducting order, so is not 
properly modelled by the random-field XY 
model. In the superconductor, this melt- 
ing transition is first-order; the XY model 
shows no such transition. However, for the 
transition from the pinned vortex lattice to 
the vortex glass, both phases have an off- 
diagonal long-range (superconducting) order, 
so the superconducting order could be merely 
a bystander, and the random-field XY model, 
which ignores this order, might capture the 
essential physics of this transition. 

There is some evidence that there may 
be two distinct superconducting phases in 
clean crystals of YBCO and BSCCO. At 
low applied magnetic fields, transport mea- 
surements show a first-order melting transi- 
tion of the vortex lattice [16,17]. At higher 
fields, the effective random pinning appears 
stronger [18], and the superconducting tran- 



sition is continuous, as expected for the melt- 
ing of a vortex glass. For BSCCO, a neutron 
scattering study saw the lattice Bragg peaks 
in the low-field regime but not in the high- 
field regime [19]. This is all consistent with 
these samples having the pinned vortex lat- 
tice phase in the low field regime and the vor- 
tex glass phase in the high field regime. Some 
signs of a transition between these two super- 
conducting phases have recently been seen in 
the nonlinear transport properties near the 
critical current in YBCO [20]. This transi- 
tion, we argue, is characterized by the prolif- 
eration of dislocations in the lattice, reducing 
it to the amorphous vortex glass. It is this 
transition that is modelled by the random- 
field XY model we have studied here. 

In conclusion, we have performed for the 
first time extensive Monte Carlo simulations 
of the d = 2 and d = 3 random field XY 
model. In both dimensions, the spacing be- 
tween static vortices grows faster than the 
pinning length obtained from the Imry-Ma 
argument as the random field amplitude goes 
to zero. In d = 3, our results are consis- 
tent with a phase transition at nonzero crit- 
ical random field into a novel topologically 
ordered (vortex-free) phase with power-law 
decay of the spin-spin correlation function. 

We thank Sue Coppersmith for discus- 
sions. The work at TRIUMF was supported 
by the NSERC of Canada. 
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FIG. 1. The fraction of plaquettes occupied 
by vortices vs. the rms random field strength. 
Triangles are d = 3, T = 1.5; squares are d = 2, 
T = 0.7. The slopes of the lines on this log-log 
plot are 4 and 16/7. These data are from simula- 
tions of single large samples (10 5 — 10 6 spins), so 
the sample-to-sample statistical errors have not 
been properly estimated. From the cases where 
we did simulate more than one sample, we ex- 
pect that the errors are roughly at the 5% level 
for the lowest fields and smaller at higher fields. 



5 10 j- 15 20 25 

FIG. 2. The magnetization correlation func- 
tion versus distance for (a) d = 3, T = 1.5, 
and (b) d = 2, T = 0.7, for the indicated 
random field strengths. The dotted lines in 
(a) are fits to simple exponentials. The error 
bars, where shown, indicate variations between 
2-3 samples. Sample-to-sample statistical errors, 
although not properly measured, appear to be 
larger, particularly at low field and large r. 
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FIG. 3. For d = 3, T = 1.5 and H c = 1.35 
the vortex density, fy, (solid symbols) and cor- 
relation length, £, (open symbols) vs. (H — H c ) 
on a log-log plot. £ is obtained from the simple 
exponential fits in Fig. 2a. 
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